THE ENERGY DENSITY IN THE PLANAR ISING MODEL 
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Abstract. We study the critical Ising model on the square lattice in bounded 
simply connected domains with + and free boundary conditions. We relate 
the energy density of the model to a discrete fermionic spinor and compute 
its scaling limit by discrete complex analysis methods. As a consequence, 
we obtain a simple exact formula for the scaling limit of the energy field one- 
point function in terms of the hyperbolic metric. This confirms the predictions 
originating in physics, but also provides a higher precision. 
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1. Introduction 

1.1. The model. The Lenz-Ising model in two dimensions is probably one of the 
most studied models for an order-disorder phase transition, exhibiting a very rich 
and interesting behavior, yet well understood both from the mathematical and 
physical viewpoints |Bax89| IMcWu73| IPal07| . 
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After Kramers and Wannicr [KrWa41_ derived the value of the critical temper- 
ature and Onsager |Ons44j analyzed the behavior of the partition function for the 
Ising model on the two-dimensional square lattice, a number of exact derivations 
were obtained by a variety of methods. Thus it is often said that the 2D Ising 
model is exactly solvable or integrable. Moreover, it has a conformally invariant 
scaling limit at criticality, which allows to use Conformal Field Theory ( CFT) or 
Schramm's SLE techniques. CFT provides predictions for quantities like the cor- 
relation functions of the spin or the energy fields, which in principle can then be 
related to SLE. 

In this paper, we obtain a rigorous, exact derivation of the one-point function 
of the energy density, matching the CFT predictions |DMS97[ ICar84[ IBuGu93j . 
We exploit the integrable structure of the 2D Ising model, but in a different way 
from the one employed in the classical literature. Our approach is rather similar to 
Kenyon's approach to the dimer model [KenOOj . 

We write the energy density in terms of discrete fermionic spinors introduced 
in [SniiOGj . These spinors solve a discrete version of a Riemann boundary value 
problem, which identifies them uniquely. In principle, this could be used to give an 
exact, albeit very complicated, formula, that one could try to simplify - a strategy 
similar to most of the earlier approaches. Instead we pass to the scaling limit, 
showing that the solution to the discrete boundary value problem approximates 
well its continuous counterpart, which can be easily written using conformal maps. 
Thus we obtain a short expression, approximating the energy density to the first 
order. Moreover, our method works in any simply connected planar domain, and 
the answer is, as expected, conformally covariant. 

Similar spinors appeared in Kadanoff and Ceva [KaCe71 and in Mercat [MerOlJ, 
but their scaling limits with boundary conditions were not discussed before [Smi06j . 

Recall that the Ising model on a graph Q is defined by a Gibbs probability 
measure on configurations of ±1 (or up/ down) spins located at the vertices: it is a 
random assignment (a x ) xeV of ±1 spins to the vertices V of Q and the probability 
of a state is proportional to its Boltzmann weight e~@ H , where f3 > is the inverse 
temperature of the model and H is the Hamiltonian, or energy, of the state a. In 
the Ising model with no external magnetic field, we have H := — <J i (T ji where 
the sum is over all the pairs of adjacent vertices of Q. 



1.2. The energy density. Let 17 be a Jordan domain and let Q$ be a discretization 
of it by a subgraph of the square grid of mesh size S > 0. We consider the Ising 



model on the graph Q$ at the critical inverse temperature j3 c — ^ In (v2 + lj; on 
the boundary of £1$, we may impose the value +1 to the spins or let them free (we 
call these + and free boundary conditions respectively) . Our main result about the 
energy density is the following: 



Theorem. Let a € £1 and for each 5 > 0, 
Then, as 6 — > 0, we have 



let (xg,ys) be the closest edge to a in 



E, 



_V2 

T xs a vs 2 



V2 



In (a) 
2ir 

In (a) 
2ir 



6 + o(S), 
S + o(S), 
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where the subscripts + and free denote the boundary conditions and Iq is the ele- 
ment of the hyperbolic metric o/fi. 

A precise version of this theorem in terms of the energy density field is given in 



Section 1.4 This result has been predicted for a long time by CFT methods (see 
[DMS97, BuGu93 for instance), notably using Cardy's celebrated mirror image 
technique |Car84| . However, CFT does not allow to determine the lattice-specific 
constant ^ appearing in front of the hyperbolic metric element. 

This is one of the first results where full conformal invariance (i.e. not only 
Mobius invariance) of a correlation function for the Ising model is actually shown. 
The proof does not appeal to the SLE machinery, although the fermionic spinor that 
we use is very similar to the one employed to prove convergence of Ising interfaces to 
SLE(3) [ChSmOQ]. A generalization of our result with mixed boundary conditions 
could also be used to deduce convergence to SLE. 

In the case of the full plane, the energy density correlations have been studied 
by Boutillier and De Tiliere, using dimer model techniques [B0DTO8, BoDT09j. 
However, their approach works in the infinite- volume limit or in periodic domains 
and does not directly apply to arbitrary bounded domains. 

In the case of the half-plane, the energy density one-point function has been 
recently obtained by Assis and McCoy [As Mcll L using transfer matrix techniques. 

The strategy for the proof of our theorem relies mainly on: 

• The introduction of a discrete fermionic spinor, which is a complex defor- 
mation of a certain partition function, and of an infinite- volume version of 
this spinor. 

• The expression of the energy density in terms of discrete fermionic spinors. 

• The proof of the convergence of the discrete spinors to continuous fermionic 
spinors, which are holomorphic functions. 

1.3. Graph notation. Let us first give some general graph notation. Let Q be a 
graph embedded in the complex plane C. 

• We denote by Vg the set of the vertices of Q, by £ g the set of its (unoriented) 
edges, by £g the set of its oriented edges. 

• We identify the vertices Vg with the corresponding points in the complex 
plane (since Q is embedded). An oriented edge is identified with the differ- 
ence of the final vertex minus the initial one. 

• Two vertices vi , v 2 € Vg are said to be adjacent if they are the endpoints 
of an edge, denoted (vi, v 2 ), and two distinct edges e\, e 2 £ Eg are said to 
be incident if they share an endvertex. 

1.3.1. Discrete domains. 

• We denote by C$ the square grid of mesh size S > 0. Its vertices and edges 
are defined by 

Vcs ■= {$ (j + ik) -j^ez}, 

£cs ■= {(vi,v 2 ) :v 1 ,v 2 G Vc s ,\vi-v 2 \ =d}. 

• In order to keep the notation as simple as possible, we will only look at finite 
induced subgraphs Q$ of C5 (two vertices of Q$ are linked by an edge in fig 
whenever they are linked in Cg), that we will also call discrete domains. 
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• For a discrete domain we denote by Q$ the dual graph of £1$: its vertices 
Vfi* are the centers of the bounded faces of Slg and two vertices of Vn* are 
linked by an edge of £q* if the corresponding faces of share an edge. 

• We denote by dVn s the set of vertices of Vc s \ Vq s that are at distance i5 
from a vertex of Vn s (i- e - that are adjacent in Cs to a vertex of Vn 5 ) and by 
d£n s C d£ c s the set of edges between a vertex of Vq s and a vertex of dVn 5 . 
The vertices in OVq s appear with multiplicity: if a vertex of Vc s \ Vq s is at 
distance 8 to several vertices of Vq 5 , then it appears as as many distinct 
elements of dVn s . In other words, there is a one-to-one correspondence 
between dVn s and d£n s - 

• We denote by dVn* the centers of the faces of C,5 that are adjacent to a face 
of Qs- We denote by d£n* C <9£cj the set of dual edges between a vertex 
of Vqj and a vertex of dVn* s . For an edge e G £n 5 we denote by e* <G £n* 5 
its dual (e and e* intersect at their midpoint). 

• We write Vn s for V f2a U <9V f2a and Vn* for V n * x U <9Vn- . 

u u u d o 

• We denote by £q C £n 5 the set of horizontal (i.e. parallel to the real axis) 
edges of fia and by £q s :— £n s \ £q the set of the vertical ones. 

• We denote by fi^ and the medial graph of Vis'- its vertices Vqm are the mid- 
points of the edges of £q s U d£a 6 and the medial edges £ n M link midpoints 
of incident edges of £ n& U d£^ s ■ 

• We say that a family (^a^o °f discrete domains (with fig C C,5 for each 
S > 0) approximates or discretizes a continuous domain if for each S > 0, 

is the largest connected induced subgraph of Cs contained in f2. 

1.3.2. Ising model with boundary conditions. The Ising model (with free boundary 
condition) on a finite graph Q (in this paper, Q will be a discrete domain Q$ or its 
dual OJ) at inverse temperature (3 > is a model whose state space is given 

by := |(c a; ) ;( . e y s : c x G {±1}|: a state assigns to every vertex x of Q & spin 

a x G {±1}- The probability of a configuration a G is 

with the energy (or Hamiltonian) Hg' ilee of a configuration cr given by 

H^ hcc (a) := - ^ crzcry, 

{x,y)eSg 

and the partition function Zg' bcc by 

Given a graph Q with boundary vertices dVg (like with <9Vn* ) the Ising model 
on Q with + boundary condition is defined as the Ising model on Q, with extra spins 
located at the vertices of dVg that are set to +1 and with energy 

H g + (<*") : = - ^2 °xO- y , 

{x,y)e£g 

where £ g is the set of edges linking vertices of Vg U dVg . 



THE ENERGY DENSITY IN THE PLANAR ISING MODEL 



5 



In this paper, we will be interested in the Ising model with free and + boundary 
conditions on discrete square grid domains Q$ at the critical inverse temperature 
f3 c := i In (v2 + l). when the mesh size 8 is small. 

We will from now on omit the inverse temperature parameter /? in the notation 
and will denote by P^ co and the probability measures of the Ising model on Q at 
j3 = /3 C with free and + boundary conditions and by E^ cc and Eg the corresponding 
expectations. 

1.4. The energy density. Let ils be a discrete domain and let as <E Vqm be the 
midpoint of a horizontal edge of fls- We introduce the two quantities (e$ (as))^ 
and (es(as))a*, called average energy density (with free and + boundary condi- 
tions), defined by 

V2 



(e s (a s ))+, := E+ 



_V2 

T n s CT ss n 



where (e$,w§) € £q s and (ng,ss) € £h* are respectively the (horizontal) edge and 
the dual (vertical) edge, the midpoint of both of which is a$ (see Figures |j~l"1and 



1.2 1. The quantity ^ is the infinite- volume limit of the product of two adjacent 
spins (it can be found in |McWu73| , Chapter VIII, Formula 4.12, for instance). The 
energy density field is the fluctuation of the product of adjacent spins around this 
limit: it measures the distribution of the energy H among the edges, as a function 
of their locations. We are considering horizontal edges on £1$ and vertical edges 
on f2J for concreteness and for making the notation simpler, but our results are 
rotationally invariant. 

We can now state the main result of this paper, which is the conformal covariance 
of the average energy density: 

Theorem 1. Let Q, be a C l simply connected domain and let a G Cl. Consider 
a family {£ls)s>o °f discrete domains approximating Q and for each S > 0, let 
a$ G Vqm be the midpoint of horizontal edge that is the closest to a. Then as 
5 — > 0, uniformly on the compact subsets of £1, we have 

I ( e «5 (as))™ ^ ( a ) ' 

6 Ztt 



s Z7T 

in (a) being the hyperbolic metric element of Q at a. Namely, in (a) := 2ijj' a (a), 
where ip a is the conformal mapping from fl to the unit disk B:={zeC:|z|<1} 
such that ip a (a) =0 and tp' a (a) > 0. 



The proof will be given in Section 1.6 



Corollary 2. The conclusions of Theorem^hold under the assumption that il is 
a Jordan domain. 

Proof. We have that (e^ (as))^ and (e,5 (a^))^ c are monotone respectively non- 
increasing and non-decreasing with respect to the discrete domain ilg, as follows 
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easily from the FKG inequality applied to the FK representation of the model (see 
| GriOG] Chapters 1 and 2, for instance): if C ilg, we have 

If O is a Jordan domain, we can approximate £1 by monotone increasing and de- 
creasing sequences of smooth domains, for which Theorem [T] applies, and deduce 
the result for ft. □ 

The central idea for proving Theorem [T]is to introduce a discrete fermionic spinor 
which is a two-point function /q s (a, z); it is defined in the next subsection. We then 
relate /n 5 to the average energy density and prove its convergence to a holomorphic 
function /q. 

1.5. Contour statistics and discrete fermionic spinor. 

1.5.1. Contour statistics. Let fl$ be a discrete domain. We denote by Cq s the set of 
edge collections uj C £n s such that every vertex v £ Vq s belongs to an even number 
of edges of uj: in other words, by Euler's theorem for walks, the edge collections 
cj £ Cfi e are the ones that consist of edges forming (not necessarily simple) closed 
contours. For an edge e £ £q s , we denote by the set of configurations w £ Cq s 

that do not contain e and by Cq 6 ^ the set of configurations that do contain e. 

Set a := v2 — 1. For a collection of edges u> C £q s , we denote by \u>\ its 
cardinality. For e £ £q s we define 

r{e+} . = \- a H 



Sis / I 



1"+} 



zr } == e « M - 



uiec 



We now have the following representation of the energy density in terms of contour 
statistics. 

Proposition 3. Let e £ £q s be a horizontal edge and let its midpoint be a £ Vqm . 
Then we have 

U Z fe +} Z n s } V% z o° +} V2 + 2 
{£s{a)) ^ ~ -7^-^-^~ 2 ^ s 2—' 

Proof. From the low-temperature expansion of the Ising model (see |Pal07j , Chapter 
1, for instance), there is a natural bijection between the configurations of spins a 
on Vfi» with + boundary condition on dVn* and the edges collections uj £ Cq 5 : one 
puts an edge e £ £n s in the edge collection uj if the spins of a at the endpoints 
of the dual edge e* £ £q* are different. It is easy to see that the probability 
measure on Cn 5 induced by this bijection gives to each edge collection uj £ Cq s 

a weight proportional to (e" 213 ) , hence to a'"', where a = \/2 — 1 as above 
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Figure 1.1. The Ising model on Vn* with + boundary condition, 
with the contours corresponding to its low-temperature expansion. 




Figure 1.2. The Ising model on Vn s with free boundary condition. 

(since /? = /3 C = iln(\/2 + l))- The event that the spins at two adjacent dual 
vertices x, y £ Vn* are the same (respectively are different) corresponds through 

the natural bijection to Cq (respectively Cq ) , where e € £q jS is such that 

e* = (x,y). Using that Zj^ + ^ + ^ = Zq 6 , we deduce the first identity. 

From the so-called high-temperature expansion (see |Pal07j . Chapter 1) we have 
that for the Ising model on Q$ with free boundary condition, the correlation of two 
spins z\ , Z2 G Vq s is equal to 

E. eCn ,(tanh(/3))l"l ' 
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where Cn s (zi,Z2) is the set of edge collections Cj such that every vertex in Vn 6 \ 
{z%, z 2 } belongs to an even number of edges of Cj and such that z%, Z2 both belong 
to an odd number of edges of Cj. At f3 = f3 c , we have tanh (/3) = a (the fact that 
tanh/3 c = e~ 2/3c actually characterizes (3 C ). 

Let us now take z±, adjacent, set e := (zi, Z2) G £n s > an d denote by Cq s (zi, Z2) 
and Cq 5 (z\, z%) the sets of Cj G Cq s {z\, Z2) such that e G Cj and e ^ Cj respectively. 
From each Cj G (z%, Z2), we can remove e and obtain an edge collection in 
(this map Cq s (21,2:2) — > C^ s + ^ is bijective) and to each Cj G Cq (2:1,22), we can 
add e and obtain an edge collection in Cq 6 ^ . Hence we have 

free Ezec+ ( Zl ,z 2 ) aH Euec n ( Zl ,z 2 ) aH V2 

Ma)) ^ = z n , + z^ V 



■'a 

^ Zo 5 2 



z n 6 2 

Using the relation Z^ + ^ + Z^f, ^ = Zj-^ , we obtain the second identity. □ 

1.5.2. Discrete fermionic spinor in bounded domain. If a G V^ju is the midpoint of 
a horizontal edge e\ G £q b and z G Vqm is the midpoint of an arbitrary edge e2 G 
£fi 5 U we denote by Cq s (a, z) the set of 7 consisting of edges of £q s \ {e 1; 62} 

and of two half-edges (half of an edge between its midpoint and one of its ends) 
such that 

• one of the half-edges has endpoints a, a + | ; 

• the other half-edge is incident to z; 

• every vertex v G Vn s belongs to an even number of edges or half-edges of 
7- 

For a configuration 7 G Cn s (a, z), we call admissible walk along 7, a sequence 
eo, ei, . . . , e„, such that 

• eo is the half-edge incident to a; 

• e„ is the half-edge incident to z; 

• ei, . . . , e„__i G £0,5 are edges; 

• ej and e^+i are incident for each j G {0, . . . , n — 1}; 

• each edge appears at most once in the walk; 

• when one follows the walk and arrives at a vertex that belongs to four edges 
or half-edges of 7 (we call this an ambiguity) , one either turns left or right 
(going straight in that case is forbidden). 

It is easy to see that, for any 7 G Cn a (a, z), such a walk always exists, though in 
general it is not unique (see Figure 1.3 1. 

Given a configuration 7 G Cq 5 (a, z) and an admissible walk along 7, we define 
the winding number W of 7, denoted W (7) G R/4-7rZ, by 

7T 

w (7) := 2 ( n < - ' 

where and n r are respectively the numbers of left turns and right turns that the 
admissible walk makes from a to z: it is the total rotation of the walk between a 
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z 



a 



Figure 1.3. An admissible walk in a configuration in Cq s (a,z). 
There are sixteen different choices of such walks in this configura- 
tion. 



and z, measured in radians. More generally, we define the winding number of a 
rectifiable curve as its total rotation from its initial point to its final point, measured 
in radians. The following lemma shows that the winding number (modulo 4tt) of a 
configuration 7 e Cq s (a, z) is actually independent of the choice of the walk on 7. 

Lemma 4. For any 7 £ Cq s (a,z), the winding number W (7) G R/47rZ is inde- 
pendent of the choice of an admissible walk along 7. 

The proof is given in Appendix A. 

Thanks to Lemma [4] we can now define the discrete fermionic spinor /n { that 
will be instrumental in our studies of the energy density. 

Definition 5. For any midpoint of a horizontal edge a € Vqm, we define the 
discrete fermionic spinor /n 5 (a, •) : V^a/ — > C by 

5 ieCn s (a,z) 



fn s (a, a) 



Z 



! n 5 



where I7I denotes the number of edges and half-edges of 7, with the half-edges 
contributing | each. 

In this way, z 1— >■ f^ s (a, z) is a function, whose value at z — a gives (up to an 
additive constant) the average energy density at a. As we will see, moving the 
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point z accross the domain will allow to gain information about the effect of the 
geometry of the domain on the energy density. 

1.5.3. Discrete fermionic spinor in the full plane. As mentioned above, the ^ 



appearing in the definition of the average energy density (Section 1.4) is the infinite 
volume limit (or full-plane) average product of two adjacent spins, which one has 
to subtract in order for the effect of the shape of the domain to be studied. We 
now introduce a full-plane version of the discrete fermionic spinor, whose definition 
a priori seems quite different from the bounded domain version. It will allow to 
represent the energy density (with the correct additive constant) in terms of the 
difference of two discrete fermionic spinors. 

Definition 6. For a,z£ Vqm with a ^ z and a being the midpoint of an horizontal 
edge, define fc s by 



f Cl (a,z) := 2cos - eT C v - 21 , - + C 



2(q-f) 2z x 
S ' s 




where Cq (zi,Z2) ■— Co (0,^2 — Zi) is the dimer coupling function of Kenyon (see 
|Ken00| ), defined by 



C(0, x + iy) := 



1 rlK p2ir e i(xe-y4>) 



4tt 2 J q J 2i sin (9) + 2 sin (<j>) 
We set / c , (a, a) := 

From our definitions up to now and from Proposition [3] we deduce the following: 

Lemma 7. Let fig be a discrete domain and a £ Vqm be the midpoint of a horizontal 
edge ofQ$. Then the average energy density can be represented as 

(e 5 (a))+, = 2(fa s -f Cs )(a,a), 

( e s( a ))nT = - 2 {fn„-fc s )(a,a). 

1.6. Convergence results and proof of Theorem [l| The core of this paper 
is the convergence of the discrete fermionic spinors to continuous ones, which are 
holomorphic functions. Let us define these functions first: for a, z € Q with a =/= z, 
we define 

fn(a,z) := ^V^M) V^¥) ^ t 1 , 

fc{a,z) := 7 r— r, 

2ir [z — a) 

where ip a is the unique conformal mapping from Q to the unit disk D with %l) a {<i) — 
and ip' a (a) > (this mapping exists by the Riemann mapping theorem). Note that 
z 1— > /n (a, z) and zi->/c (a, z) both have a simple pole of residue j- at z = a and 
that (fa — fc) (a, z) hence extends holomorphically to z — a. 
We can now state the key theorem of this paper: 
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Theorem 8. For each S > 0, identify a G C with the closest midpoint of a hori- 
zontal edge of C$ and z €E C wii/i the closest midpoint of an edge ofC$. Then we 
have the following convergence results 



fa s (a, z) t-> fn (a, z) Va, z € £l : a ^ z, 



1 

xfc s (a,z) — » fc{a,z) Va, zeC:a^z, 

O o->0 

7 (/n 6 - /cj (a, z) (/n - fc) (a, z) Va, z € fi. 

S-)-Q 

where the convergence of \fn s is uniform on the compact subsets of ft x il away 
from the diagonal {(w, w) : w £ fi}, the convergence of \fc e * s uniform on C x C 
away /rom £/ie diagonal {(w,w) : to € C} and the convergence of \ (fn 5 ~ fcs) * s 
uniform on the compact subsets ofQxtt. 

From this result, the proof of the main theorem follows readily: since we have 
(Lemma [7]) 

(e<5 (a)) n* = 2 (/n 5 - /c s ) 0, a) , 

( e <5 ( a ))o° C = - 2 (fas - fas) (a, «) ) 
and since | (/n 4 — fc s ) converges to (fn — fc) (a, a), it suffices to check that 

V^W^M ^I 1 - — ) -> V£ (a) , 

(Zj Z - a/ z^a 

which follows readily by verifying that 



^U^)V^Azj i ^ 

( z ) Z — a z^a 

To prove this, notice that since 

— ; — = h A + CJ (z — a) z — » a, 

Va (z) 2 - a 

by squaring this expression, it suffices to check that the residue of ip' a (z) /ip^ (z) 
at z = a vanishes. By contour integrating on a small circle C around a and using 
change of variable formula (since ip a is conformal), we have 

Ani - A £ Vo ( z ) dz £ dw 



1.7. Proof of convergence. The rest of this paper is devoted to the proof of the 
key theorem (Theorem j8|. This proof consists mainly of two parts: 

• The analysis of the discrete fermionic spinors fa s and fc 6 as functions of 
their second variable (Section 2): 

— We prove that fa s (a, •) and fe s (a, •) are discrete holomorphic (in a 
specific sense) on Vqm \ {a} and V^m \ {a} respectively (Propositions 



13|and|14|. 

— We show that /q 5 (a, •) and fc e (a, ■) have the same discrete singularity 
at a: their difference (fn s — fc s ) ( a i ') is hence discrete holomorphic on 



Vqm (Propositions 15 16 and 17) 
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— We observe that /q 4 (a, •) has some specific boundary values on OqVqm 
(Proposition 18 1. 

The proof of convergence of functions \fa e , \fc s , \ (fn s ~ fc s ) (Section 
3): 

— The convergence of \fc 5 follows directly from the convergence result 
of Kenyon for the dimer coupling function (Theorem 25 I . 

— We show that the family of functions (| (/n 5 — fcs ))s>o * s P recom P ac t 
on the compact subsets of £ 1 x Q: it admits convergent subsequcntial 
limits as 8 — > (Proposition 26 1. Hence (^fn s ) s>0 is also precompact 
on the compact subsets of f2 x Q away from the diagonal. 

— We identify the 5 — > subsequcntial limits of j/n 5 with the function 
/n (Proposition 30 1 . This allows to conclude the proof of Theorem M 
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2. Analysis of the Discrete Fermionic Spinors 

In this section, we study the properties of the discrete fermionic spinors fo s and 
fc s that follow from their constructions. In the next section, we will use these 
properties to prove Theorem [8] 

We study both spinors as functions /q s (a, •) and fc s (a, •) of their second vari- 
able, keeping fixed the medial vertex a 6 Vqm (which is the midpoint of a horizontal 
edge of fia). 

Let us first introduce the discrete versions of the differential operators d and A 
that will be useful in this paper: for a C- valued function /, we define, wherever it 
makes sense (i.e. for vertices of Vn 5 U Vnj), 

WW - /H)-/H)+.(/H)-/H)). 

A 4 / (x) := f(x + 6) + f(x + iS) + f(x-S) + f(x- i6) - 4/ (x) . 

In the case where one has that a vertex y € {x ± S, x ± id} belongs to dVn 5 in the 
definition of A<5, the boundary vertex is the one identified with the edge (x, y) £ 
dS Qs . 

If e = xy € £n 5 IS an oriented edge with x, y € Vn s and / is a function Vq 5 — > C 
we denote by d e f the discrete partial derivative defined by d e f := / (y) — f (x). 

2.1. Discrete holomorphicity. It turns out that the functions fa s (a, ■) and fc s (a, ■) 

are discrete holomorphic in a specific sense, which we call s-holomorphicity or spin- 
holomorphicity. 
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FIGURE 2.1. The lines associated with the medial edges of 

Let us first define this notion. With any medial edge e € £ c m , we associate a 
line £ (e) C C of the complex plane defined by 

£(e) := (d-u)~»]R= {(d - v)~h : t e r| , 

where u g Vc s is the closest vertex to e and d € Vc* is the closest dual vertex to e. 
On the square lattice, the four possible lines that we obtain are e ± ^IR and e ±3 ff 1 ]R. 
When I := e ,e ]R is a line in the complex plane passing through the origin, let us 
denote by Pe the orthogonal projection on £, defined by 

P e [z] := - (z + e 2i9 z) Vz G C. 

Definition 9. Let Ais C V c m be a collection of medial vertices. We say that 
/ : Ais — > C is s-holomorphic on Ms if for any two medial vertices x, y € .A/Ja that 
are adjacent in C^ 1 , 

P, (e) [f{x)] = P t[e) lf(y)}, 

where e = (x,y) £ S^m. 

Remark 10. Our definition is the same as the one introduced in [SmilO , except 
that the lines that we consider are rotated by a phase of eT, and that our lattice 
is rotated by an angle of ?, compared to the definitions there. In [ChSm09 this 
definition is also used, in the more general context of isoradial graphs. 

Remark 11. The definition of s-holomorphicity implies that a discrete version of 
the Cauchy-Riemann equations is satisfied: if / : Ais C V c a/ — > C is s-holomorphic 

and v £ Vc 5 U Vcj is such that the four medial vertices v ± |,w ± y are in Ais, 
then we have: 

d s f(v)=0. 

This can be found in [SmilO] (it follows by taking a linear combination of the 
four s-holomorphicity relations between the values / (y ± f ) , / (y ± y)), as well 
as the fact that satisfying this difference equation is strictly weaker than being 
s-holomorphic. 

Remark 12. If Xs — |v,Wi+i € £qm : i G Z/nz| is a simple counterclockwise-oriented 
closed discrete contour of medial edges and As is the collection of points in Vn s UVn» 
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surrounded by \$, then it is easy to check that for any function / : V&,m — > 
have 



we 



E 



f(Vi) + f(v i+1 ) 



(v i+ i - Vi) = iS ^ d $f ( v ) 

z6A s 



Viv i+1 e\s 

In particular this sum vanishes if / is discrete holomorphic. 

Proposition 13. The function /q s (a, •) is s-holomorphic on Vqm \ {a}. 

Proof. Let z,w G Vqm \ {a} be two adjacent medial vertices and let e G Sqm be 
the medial edge linking them. Suppose that z is the midpoint of a horizontal edge 
and that w is the midpoint of a vertical edge. We prove the result in the case where 



w = z + 



l+i 



S (the other ones are symmetric) . Denote by ft the half-edge between z 
and z+ 1 G Vn 5 and by ft the half-edge between z+f and w. For any 7 g Cq s (a, z). 
define (7) := 7 © ft © ft as the symmetric difference of 7 with if ft is not 

in 7, add it, otherwise remove it, and similarly for ft. Clearly, (p is an involution 
mapping Cq, s (a,z) to Cn s (a,w) and vice versa. Moreover, for 7 € Cn« (^j^)j we 
have 



(2.1) 



|w( 7 ) 



a lv(7)l e -fW( V (7)) 



This identity follows from considering the four possible ('MSGS j ctS shown in Figure 
I2~2l 

(1) If ft £ 7 and ft ^ 7, then we have e ~2 w(7) € K, (t) I = M + 1 and 

e -*W(^( 7 )) = e -^ e -fW( 7 ), 

(2) If ft ^ 7 and ft G 7, we have e -f w (7) g |^( 7 )| = | 7 | and e -f w M7)) = 
e~^" e~%~ vv ( , y) : there are a number of subcases, as shown in Figure 2.2 for 
which these relations are satisfied. 

(3) If ft G 7 and ft G 7, we have e~i w W G iR, \<p(i)\ = \j\ - 1 and 
e-| w (cW) = e -Te~5 w h); in this case, we can always choose an ad- 
missible walk on 7 that is like in Figure |2~2| 

(4) If ft G 7 and ft £ 7, we have e"5 w (T) G iR, |<£ (7)! = M and e"i w ^(t)) = 

In all the four cases, it is then straightforward to check that equation |2.1| is satisfied. 
By definition of /n s (Section [5|, we finally deduce 

1 



P.-affi.t/n, (o,z)] 



1 



E 



a l7l e -|W( 7 ) 



a l7l e -fW( 7 ) 



which is the s-holomorphicity equation. □ 

The full-plane discrete spinor is also s-holomorphic: 

Proposition 14. The function fc s (a, ■) is s-holomorphic on V c m \ {a}. 

Proof. This follows directly from the definition of fc s (Definition |6| and of Lemma 
[32] in Appendix B. □ 
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Figure 2.2. The four possible cases in the proof of Proposition 
1 1 3 1 how a configuration is changed after two half-edges between z 
and w are removed or added. 
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2.2. Singularity. Near the medial vertex a, the functions fa s (a, •) and fa s (a, •) 
are not s-holomorphic: they both have a discrete singularity, but of the same nature, 
and consequently the difference (fa s — fa 5 ) (a, ■) is s-holomorphic on Vn m , including 
the point a. Rather than defining the notion of discrete singularity let us simply 
describe the relations that these functions satisfy near a. For x E {±1 ± i}, let us 



denote by a x :- 



x5 



£ Vqm the medial vertex adjacent to a, by e x & Sqm the 



medial edge between a and a x and by £ x the line £(e x ) (as in Definition |9j) . Let 
e £ £q 5 be the horizontal edge with midpoint a. Recall that fa s (a, a) is defined as 

Z^ g +} /Zn s (see Definition [i]). 

Proposition 15. Near a, the function fa s (a, •) satisfies the relations 

P^i +i [fa s (a, a)] = P^i+i [fa s (a, a 1+i )] , 
[fa s (a, a)} = Pe^ [fa s (a, ai_i)] , 



Pi- 1+t [fa s (a, a) - 1] 

Pl- X -i [fas ( a : a ) - 1] 



Pt- 1+l [fas (a,a_i+i)] 
P^_!_i [/n s (a,a_i_i)] 



Proof. The first two relations are the s-holomorphicity relations and they are ob- 
tained in exactly the same way as the s-holomorphicity relations away from a. 
Indeed, let us take the same notation as in the proof of Proposition 13 and consider 
the involutions (p 1+i : C^ +} -t C ns (o,a 1+i ), and ipi_i : C^ +} -t Cn 5 (a, ax-i) de- 
fined by <p 1+i (7) := 7© (a,ax) © (ax,ax+i) and </?i_ 4 (7) := 7© (a, ax) © (ai,ai_,) 
respectively: as in the proof of Proposition |15| we have that these involutions pre- 
serve the projections on £x+i and lx~i respectively (a configuration 7 e IS 
interpreted as a configuration of winding 0). 

For the last two relations, we have that the involutions tp—x+i '■ 



and : CX -> C f2s (a,a_i_ i ), respectively defined by <y9_ 1+i (7) := 7 



(o,a_i) © (a_i,a_i +i ) and y>-i-i (7) := 7 ' 
for any 7 € C^ e , we have 



(o,o_i) © ( 







= P^ 1+l 




iw(-r) 1 
2 




>l" 






iw(-r) i 
2 



<x_i, a_i_i 



a l»»-l+i(7)l e - 



Cn 5 (a, a_x 

: 7 © 
are such that 



iW(|p_ Hi h)) 



lv-i-i(7)l, 



iw( y _ 1 _ i ( T )) 



where 7 is interpreted as a configuration with a path from a to a that makes a loop, 
of winding nu mbe r ±2ir. This follows from the same considerations as in the proof 
of Proposition 13 Hence, since ^ /Zq 5 = 1 — /n s ( a i a ) ; we obtain, by summing 



le— 1 

the above equations over all 7 £ Cq s , the last two identities of Proposition 15 □ 

The function fa s (a, ■) has the same type of discrete singularity as fa 5 (a, ■): 

Proposition 16. Near a, the function fa s (a, •) satisfies exactly the same projection 
relations as the ones satisfied by the function fa s (a, •), given by Proposition 15 



□ 



Proof. See Appendix B, Proposition |33| 

From Propositions [T3| [14| [15] and [T6| we readily deduce the following: 
Proposition 17. The function (fa s ~ fas) ( a r) : Vq m — > C is s-holomorphic on 
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Figure 2.3. When z G OqVqm , the winding number of the walks 
on Cq s (a, z) is determined modulo 2tt. 



2.3. Boundary values. A crucial piece of information to understand the effect of 
the geometry of the discrete domain VL$ on the average energy density at a G Vqm 
is the boundary behavior of /q s (a,-). On the set of boundary medial vertices 
8oVqm, which link a vertex of tt$ and a vertex of dQ$, the argument of /n 5 (a, •) is 
determined modulo n. For each z G OoVqm , with z being the midpoint of an edge 

e G 9£q s between a vertex x G Vq s and a vertex y G dVj^ , denote by ^ out (z) G d£^i s 
the oriented outward-pointing edge at z, identified with the number y ~ x: it is a 
discrete analogue of the outward-pointing normal to the domain. 

Proposition 18. On dgV^M , the argument of the value of /q s (a,-) is determined 
(modulo ix): for each z G OqVqm , we have 

3m (/n s (a, z) i/ D 3 ut (z)^) = 0. 
Proof. From topological considerations, we have that if z G OqVqm and 7 G Cq 6 (a, z). 
then 3m (e~i w ^ v^ ut (z)\ = 0: the winding number of any admissible walk from 



a to z is determined modulo 2-k (see Figure 2.3 1 and it is easy to check that e 2 w (t) 



is a real multiple of {v OVL t ( z )) 2 ■ Hence, the result follows from the definition of 
fas- □ 

Remark 19. This is the same kind of Riemann-type boundary conditions as in 
|SmilO|ICTiSm09"] . Notice that in these papers, the argument of the function on the 
boundary is fully determined (not only modulo 7r). 

2.4. Discrete integration. An essential tool that we will use for deriving the 
convergence of fa s (a, •) is the possibility to define a discrete version of the anti- 
derivative of the square of an s-holomorphic function, cf. |Smil0j . 

Proposition 20. Let f : V p m — > C be an s-holomorphic function on a discrete 
domainT>$ and letx G V-p s UV-p* (where Vx> 6 = VT> s UdV-D 5 andVv* — Vx>*U9Vp* ). 
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Then there exists a (possibly multivalued) discrete analogue l Xj g [/] : Vx> s UVx>j 
of the antiderivative 

-*(// 

uniquely defined by I Xt g [/] (x) := and 

l x , S [f](b)-l x , s [f](w) = V^SlP^lfiy)}] 2 = V26\P e{e * ) [f(z)}\ 2 

V6 G Va { , w G Vnt : \b-w\ = —, 

where e* — (y,z) G £p,M * s ^ e medial edge which is between b and w. If T>s is 
simply connected, then the function l Xt s [/] is globally well-defined (single-valued) . 
When the choice of the point x is irrelevant, we will ommit it and simply write 
hlf}- 

Remark 21. It follows from the definition of 1$ [/] that for any pair of adjacent 
vertices x, y G Vv s , we have 

I, [/] (x) h [/] (y) = -«e (f 2 (y - x)\ , 

and similarly if x, y G Vpj are adjacent dual vertices. From there it is easy to see 
that if the mesh size is small, lg [/] is a good approximation of — 3?e (j / 2 ). 

We denote by I* s [f] and 1° s [/] the restrictions of I x ,s [f] to Vp a and Vx>* 
respectively. We have the following: 

Proposition 22. The function I* [/] : Vz> 5 — > K is discrete subharmonic and the 
function \° s [/] : V-p* — > K is discrete superharmonic: we have 

a 4 i; [/](«) > o vueVp,, 

Aal? [/](«) < VuGVpj. 
if m G dV-pM , then we have 

^out(m) 1 * [/] = Sm (/ M V*ut M) ~ (/ M ^out > 

where z/ out (m) G 9£x> 5 * s ^ e oriented edge from a G 2?^ to b £ dT>s, the midpoint 
of which is m. 

Proof. For the subharmonicity/superharmonicity deduced from the s-holomorphicity 
of /, see Lemma 3.8 in |SmilO| (the fact that the phases are different does not affect 
the result) . The normal derivative statement follows directly from the definition of 
h [/]. □ 

In the case of the discrete fermionic spinor /q s (•, •), the boundary condition for 
I<5 [fus ( a i ')] (') becomes particularly simple. 

Proposition 23. The function I| [fn s (a, •)] : Vq* — > K is constant on <9V$i* and 
for each m G OqV^ 



6 



^ ont (m)lj [fn e (a, •)] = - \h 5 {a,m)Y 
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Proof. The first statement follows from the construction of Ig [fn e (a, •)] and from 



the boundary condition for /q 4 (Proposition 18). 

The statement for I* [fn s (a, •)] follows directly from Proposition 22 and the 



boundary condition for f^ 5 (Proposition 18 again). □ 



Remark 24. Note that lg [fn e i a , •)] is single-valued (as follows from Proposition 23) 
and well-defined on Vn s U Vn* but that the presence of a singularity near a implies 
that and K [fn 6 (a, ■)] and I| [/q 4 (a, •)] are (at least a priori) not subharmonic or 
superharmonic near a (more precisely at a ± | , a ± ^ ) . 

3. Convergence of the Discrete Fermionic Spinors 

We now turn to the convergence of the three functions j/n a , |/c 5) j (/n a — /c s ) 
as <5 — > (Theorem |8|. For this, we use the discrete results derived in the previous 
section: the s-holomorphicity, the discrete singularity and the boundary values. As 
we will discuss convergence questions, we will always, when necessary, identify the 
points of the complex plane with the closest vertices on the graphs considered. In 
this way, we will extend functions defined on the vertices of the graphs fls, fig, fig 1 
to functions defined on fl. In particular, for the discrete holomorphic spinors, when 
we write fn s (a, z) or fc s (a, z) for a, z £ Q, we identify a with the closest midpoint 
of a horizontal edge of Sq and z with the closest midpoint of an arbitrary edge of 

The convergence of fe 6 almost immediately follows from the work of Kenyon 
[KenOOj : 



fc s (a, z) T~t fc ( a ' z ) 



Theorem 25. For any e > 0. we have 

1 

uniformly on {(a, z) £ C 2 : \a — z\ > e}, where 

fc(a,z) = — — r. 

Zir [z — a) 

Proof. See the last paragraph of Appendix B. □ 

For the convergence of \fn 5 and I {fu 5 — fc s )> we proceed in two steps: we 
first show that the family of functions Q (fn s — fc s ))g >0 i s precompact. Precom- 
pactness for \^fns)g > o w ^ then readily follow from Theorem 25 We then iden- 
tify uniquely the subsequential limits of (^fn 5 )g >0 ', this also identifies the ones of 
{g (fcis ~ fcs)) g > Q- 

3.1. Precompactness. We now state our main precompactness result: 
Proposition 26. The family of functions 

6 Js>o 
is precompact in the topology of uniform convergence on the compact subsets of 
fl x fl, and hence the family of functions 



(a,z) ^ -fa e (a,z) 



s>o 
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is precompact in the topology of uniform convergence on the compact subsets of 
Q x Q that are away from the diagonal. 



Proof. Set f n s := /n 6 — fc s - By Proposition 18 we have that for any x £ 8qVqm, 



(/§ ( a > x ) ■ "out (»)) = (/c« (a, a:) • ^ 2 ut (x 



By Theorem 25 and the fact that dfl is smooth (and hence the number of medial 
vertices in <9 Vqm is O (<5 -1 )) we deduce that the family of functions 



a i — y 



\ 



5>0 



is uniformly bounded and equicontinuous on the compact subsets of fi (since y/c a 
is uniformly convergent by Theorem 25). By Proposition 27 below, we obtain that 
the family of functions 



JQs 



v nf x Vqm 



5>0 



is uniformly bounded and equicontinuous and hence we get the desired result by 
extending the functions /q s in a uniformly continuous way to f2 x fl (for instance 
by piecewise-linear interpolation) and by using then Arzela-Ascoli theorem. □ 

Proposition 27. There exists a universal constant C > such that for each 5 > 
and any s-holomorphic function us : Vqm — ¥ C, we have, for any v £ Vqm \9oVj 2 a/ , 



M«)l < c 



\ 



1 



\\Vsu s (v) 



C 





(u s (x) ■ v* ut 


(x) 


\ 2 

) 6 


dist | 








6 


(u 5 (x) ■ V^ut 


(x) 


) 6 



^ diat(v,d V Q Mj 
where (v ) = (wa (u + 6) — us (v) ,us (v + iS) — us (v)). 

Proof. Consider the function Is [us], normalized to be at an arbitrary point. By 
subharmonicity (Proposition 22 1, a discrete integration by parts, and again by 
Proposition |22| we obtain 

< Yl A sl' s [us}= Y ^o„ t (,)I*M 



Y ( 5m ( US ( X ) ' "out ( x ) 



(Kc (us (x) ■ ^ G 2 ut (x) 



and from the last identity we also deduce that 

| 2 <2 Y, (%rti(us(x)-vi t (x) 



Y \ us ^ 



x<£d V oM 
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On the other hand, from the construction of Is [u], it is easy to see that 

\ u s(x)\ 2 -S 

By subharmonicity of I* [us], superharmonicity of I| [us] and the construction of 



max \l s [u s ] (z)\ < 

zedVn 6 UdV n * 



Is [us] (Proposition 20 ) , we have 



niax_ |I 5 [us] (w)\ = max \I S [u s ] (z) | . 

weVn s UV Q * z£dVa s UdV n * 

By Theorem 3.12 in |ChSm09] (the construction there is the same as the one of 
our paper, up to a multiplication by an overall complex factor, which does not 
affect the result), there exists then a universal constant C > such that for any 
v E Vqm \ d V Q M 

. max «,eVo,uv n . I^MHI 

MtOf < C ^ 

dist v, OoVqm 



max 



V^uVo. I 1 * M HI 



weVn 

||V^»|| 2 < C '- 

dist ( v, d V n * 

We therefore deduce the desired inequalities. □ 

3.2. Identification of the limit. We can now uniquely identify the subsequential 
limits of (j/n 5 ) (5>0 as 5 — > (we will often make a slight of abuse of notation and 
simply denote the family of functions by j/n a )- Let us start with a characterization 



of the continuous fermionic spinor /q (a, •) (defined in Section 1.6). 
Lemma 28. The function /q (a, •) is the unique holomorphic function such that 

z ^ f n (a, z) - -—— 

27r [z — a) 

is bounded near z = a and such that 

(3.1) 9m (/n (a, z) vL («)) = Vz e dQ, 

where v> out denotes the outward-pointing normal to dQ. 

The boundary condition \3.1\ is equivalent to the condition that the antiderivative 
F (z) — —die ( J z /q (a, w) dw) is single-valued on Q \ {a}, constant on 9f2 and 
satisfies 

d Vont{z) F < o Vz e dn, 

where d Vovit ^F denotes the normal derivative of F in the outward-pointing direc- 
tion. 



Proof. It is straightforward to check from the definition (Section 1.6| that /n (a, •) 
has a simple pole of order 1 and residue at z — a and satisfies the boundary 
condition 13.11 

Let / be another function with the same pole and boundary condition. Then 
the function g defined by g (z) :— /q (a, z) — f (z) extends holomorphically to f2 
and satisfies 



(p 0) ^out (- 



= Vze dCl. 
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The function G : Q — > C defined by G (w) := f g 2 (z) dz has constant real part on 
dfl and hence is constant on CI. by the maximum principle and the Cauchy-Riemann 
equations. Hence fn (a, ■) — f (•). 

For the second part of the statement, notice that the boundary condition |3.1| 
implies that /q (a, •) 2/ ou t (•) is purely real on dfl, and hence F must be constant 
on dQ. (when going along the boundary, one integrates (a, •) dr, where r is the 
tangent to the boundary, which is orthogonal to the normal ^ ou t); this implies that 
F is single-valued on \ {a} (since Q is simply connected) and that 

d„ out{z) F = -\f n (a,z)\ 2 vzeaa 

Conversely, it is easy to check that if F is constant on dft, then for any z £ d£l, we 
have (a, z) v out (z) € R. Moreover, for any z £ <9f2, we have that dv OVLt (z)F < 

implies /q (a, z) f out (z) > 0, which is equivalent to 3m ^/n (a, z) i/ Q 2 ut (z)^ =0. □ 

Let us also give a lemma which will be useful to connect the discrete spinors to 
the continuous ones: 

Lemma 29. We have the uniform bound 

l/n,(*)| 2 -*| x. 



Proof. This follows directly from the the proof of precompactness of (| (fa s — /c a )) , 




(Proposition 26 1, the convergence of ?fc t (Theorem 25 1 and the fact that is 
smooth (the number of medial vertices in 9oVn™ is O (i^ 1 )). □ 

We now pass to the identification of the subsequential limits of j/n a (<v) ^ 
Proposition 30. Let 5 n be a sequence with 5 n — > such that j^fn s (a, ■) — > 

n— >oo n n n— ¥oo 

/(•), uniformly on the compact subsets ofQ \ {a}. Then /(•) = /si (a, •), where 

/n (a, •) is defined in Section \l.6\ 

Proof. For each <5 > 0, set / 5 (•) := i/n 4 (a, •). 

Let us first remark that / (•) is holomorphic on \ {a}, as it satisfies Morera's 
condition: the integral of / (•) on any contractible contour vanishes, since it can be 
approximated by a Riemann sum involving fg n (•) (for S n small), which vanishes 
identically as explained in Remark 12 Fix a point p £ 51 \ {a}. By Lemma 
28 to identify / with fn(a,-), it suffices to check that F, defined by F (z) :— 

— 3f?e ( J* f 2 (w) du?j , satisfies the conditions of the second part of that lemma. 
Let Fs : Vn s U Vu* — > K be the discrete antiderivative I Pt s [fs], as defined in 



Proposition 20 and let F* and F$ denote the restrictions of Fg to Vsi 5 and Vn* , 



which are discrete subharmonic and superharmonic respectively (away from a), by 



Proposition 22 By Proposition 23 the function F' is constant on 8Vq* ; denote by 



Fg (Oil) this value. Fix a smootndoubly connected domain Tc!l \ {a} such that 
dil C <9T and dist (a, <9T) > (one of the component of <9Y is dil and the other 
is a simple loop surrounding a). Let us write Fg =: i/J + S* and FJ 1 =: Hg + Sg, 
where 

• H' : Vy s — > E is discrete harmonic, with H' := F 1 * on 9Vt 5 , 
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• S' : Vy 6 — > K is discrete subharmonic, with S' := on dVr s , 

• tig : Vtj — > K is discrete harmonic, with Hg := FJ" on <9Vyj, 

• S| : Vtj —> K is discrete superharmonic, with 5| := on dV-r*. 
Let us further decompose H ' as A* + B' , where 

• A° is discrete harmonic with A* := Fg (dtt) on <9Vn a and A* := H' on 
dVr 5 \dVn s - 

• £?* is discrete harmonic with B' := if* — Fg (dQ) on 9Vn 5 and B' := on 
dV Ts \dV ns . 

The situation is hence the following: for any z € Vn 4 and «; € Vn* such that 
|z — 10 1 = S/y2, from the construction of i 7 ^, the superharmonicity of Fg and the 
subharmonicity of F* , we have 

(3.2) H° s (to) < F 4 ° ( W ) < F; (z) < tf| (z) = (z) + B' s (z) . 



It follows easily from Remark 21 that as n — !• 00, we have that i 7 ^ — > F, 



2<S 



H| n and tf^ are 



and F* there, and 



uniformly on the compact subsets of Q \ {a} (since fg n —> f) . 

Let us now check that F satisfies the conditions of Lemma 
uniformly close to each other on dT \ dft (they are equal to F ; 
these functions are uniformly close to each other near dT \ <9f2, as follows easily 
from the convergence of fg n ) . To control Bg n , we use the following lemma, which 
is proven at the end of the section. 

Lemma 31. As n — > 00, B* uniformly on the compact subsets o/T. 

Observe that F$ n (dfl) is uniformly bounded. Suppose indeed that it would not 
be the case and (by extracting a subsequence) that Fg (dft) — > 00 (say). We would 
have Hg ~ > 00, since H$ is harmonic and bounded on dT \ dft (since it is equal 
to Fg there). We also would have A* — > 00, for the same reasons. By Equation 



3.2 and Lemma 31 it would imply that F$ n would blow up on T, which would 
contradict the fact that it converges uniformly to F on the compact subsets of T. 

We deduce that Hg and A' are uniformly bounded on T. 

We have that Hg —> F and A' — > F as n — >• 00, uniformly on the compact 
subsets of T. From the discrete Beurling estimate (see |Kes87j) and the uniform 
boundedness of Hg and A' near d£l, we readily obtain 



lim sup I A' Sn (z) - Fg n (dSl) \ — > 0, 

n— hoo z — 

lim sup I H° Sn (z) - Fg n (dft) I — > 0, 

and we deduce that F continuously extends to dtt and is constant there. 

To show that d Vmt ( z )F < for all z e c?Sl, we consider the harmonic conjugate C: 
it is the unique function C (defined on the universal cover of f2\ {a} and normalized 
to be at an arbitrary interior point x) such that F + iC is holormophic. By the 
Cauchy-Riemann equations, we have 

W = ^« C VzedO, 
where 9 Tccw ( z ) is the tangential derivative on dfl in counterclockwise direction and 
the condition d Vau% i z \F < becomes <9 rccw( - z )C < 0. This latter condition is equiv- 
alent to the one that C is non-increasing when going counterclockwise along (the 
universal cover of) dfl. 
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Let us now check that this condition is satisfied. Take T as before and denote 
by T its universal cover. 

For each S > 0, let C| : V-j-„ — > R be the discrete harmonic conjugate of H° 
(lifted to Vf ), defined by integrating the discrete Cauchy-Riemann equations 

d s (H' s + iC° s ) (z) = VzeV f M, 

and with the normalization C| (x) = 0. By subharmonicity of F* , we have F' < H° 
on V x m and hence, since F' = H* on 8Vqm . 

d VaM(z) H' 5 < d VaM[z) F' & < Vz e 0V nf , 

and we deduce by the discrete Cauchy-Riemann equations that C| is non-increasing 
when going along the universal cover of dV^ in counterclockwise direction. 

On the compact subsets of T, since the (normalized) discrete derivatives of H* 
converge uniformly (see Remark 3.2 in [ChSm08 ) to the derivatives of F, it is 
easy to check that C| also converges uniformly to C . Since C| is non-increasing 
(when going along the universal cover of <9f2) , we have that Cg is locally uniformly 
bounded (uniformly with respect to n) on the universal cover of dfl (if it would 
blow up there as n —> oo, it would also blow up on T), and hence it is bounded 
everywhere on the closure of T. 

From there, we deduce that C is non-increasing on the (counterclockwise-oriented) 
universal cover of 9f2: if it would not be the case, using again the discrete Beurling 
estimate [Kes87 , we would obtain a contradiction (in the n — > oo limit) with the 
fact that Cg is non-decreasing. □ 



Proof of Lemma \Sl\ For z G dVr s - let us write Pg (z, ■) : Vx^ — > K for the dis- 
crete harmonic function such that Pg (z, •) = 1{ Z } (•) on 8V^ S (this is the discrete 
harmonic measure of {z}). By uniqueness of the solution to the discrete Dirichlet 
problem, we can write 

B 5 (y)= Y, B s (z)P(z,y) Vj/eV Ti . 

As 6 — > 0, we have that Ps (x, ■) — > on the compact subsets of T, uniformly 
with respect to x (this follows directly from Proposition 2.11 in [ChSm08 ). By the 
construction of Fg and the boundary conditions (Propositions [18] and 23 1 , we have 

B S (z) = F s (z) - F s (dtl) = ^2 cos (j^J \f s (m)| 2 6 

for any z £ dVn s , where m € 8qVqm is the midpoint of the edge between z and its 
neighbor in Vn a ■ Since 

E l/*»M| a *n- 
med V n M 



is uniformly bounded by Lemma 29 we readily deduce that Bg n — > uniformly on 
the compact subsets of T. □ 
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Figure 3.1. The oriented loop L formed by adding the curve [i 
(dotted) to A. In this case, Ai = 6 and N2 = 4. 



Appendix A 



We give here the proof of Lemma [4] for a configuration 7 G Cn 5 (a, z), the 



winding (modulo 47r) of an admissible walk on 7 (see Figure 1.3 1 is independent of 
the choice of that walk. 

Proof of Lemma [7J Without loss of generality, half-edges of 7 emanate from z and 
a in the same direction, so the winding is a multiple of 2tt. 

Add a curve \x from z to a, which emanates in opposite direction from 7 and run 
slightly off the lattice, so that /x is transversal to 7 when an intersection occurs (see 



Figure 3.1 1. Let N\ be the number of intersections of /1 with 7. 

Take any admissible walk A along 7. The rest of 7 can be split into disjoint 
cycles. So, if N2 is the number of intersections of A with 7, then N2 = N-y (mod 2). 
Indeed, their difference comes from cycles, which are disjoint from A and so intersect 
fi an even number of times (see Figure |3.1[ ) . 

The concatenation of A and /i (when oriented) forms a loop C, which has sev- 
eral intersections (when A and /j, run transversally) . At each of those, change the 
connection so that there is no intersection, but instead two turns - one left and 
one right. Each of A2 rearrangements either adds or removes one loop, so after the 



procedure, C splits into N 3 simple loops with N 3 = N 2 (mod 2) (see Figure 3.2 1. 

Each of the A3 simple loops has winding 2tt or — 2ir, so W (C) = 2ttN^ = 
2nNi (mod 4tt). 

We conclude that, mod 4tt, W(A) = W(£) - W(/i) = Nx - W(/z) and so 
W(A) (mod 47r) is independent of its particular choice. □ 
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Figure 3.2. The five simple loops obtained from C after four re- 
arrangments (and discarding the loops that were not part of A). 



Appendix B 

We prove here som e tech nical results concerning the discrete full-plane spinor, 



1.5.3 



Let us denote by Cg (•, •) := C'o the rescaled 



introduced in Section 

version of Kenyon's coupling function (defined in [Ken OO ) . 
Lemma 32. With the notation and assumptions of Proposition [L/j the functions 
Gi : V c m \{a} C 

z h-> eT fc s ( a +L z ) + c s ( a ~% z 



G 2 : V c¥ \ {a} 



-> < >■ | f ( <*--*) + C s (a+~,z 



are s-holomorphic. 



Proof. Set r\ := eT and for any vertex z and any /i <E {±1, ±i}, set z^ := z + ^p. 
By translation invariance, we have 

Gi (z) = TJ (C s (a, z-i) + C s (a, Zi)) , 
G 2 (z) = ir](C$(a,Zi) + C s (a 1 z^ i )) , 

where, on the right hand sides, the two values of Cs (a, ■) are orthogonal: one is 
purely real and the other purely imaginary. Let x, y G ~^c^\{a} be two adjacent 
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medial vertices, with x being the midpoint of a horizontal edge of £q s and y the mid- 
point of a vertical one, and let e := (x,y) € £ c m . Then, there are four possibilities 
for the line I :— £ (e): 

• If £ = 77M, we have that x = y + ^p<5 and 

P e [GUx)] = r)C s (a,X- 1 )= V C 5 {a,y i ) = P e [G 1 (y)}, 
Pe [G 2 0)] = ivCs (a, X-i) = ii]C s (a, y x ) = P e [G 2 (y)} ■ 

• If £ = ry 3 IR, we have that x = y — and 

P/[Gi(ar)] = r l C s (a,x i )=r)C s (a,y- 1 ) = Pt[G 1 (y)}, 
Pi [G 2 (x)] = ii]C s (a, xi) = ii]C s (a, = P e [G 2 {y)\ ■ 

• If £ — TjM., we have x — y + ^-^S and 

Pi [Gx (x) - Gi (y)\ = —^(C s (a 1 ,x_ 1 )+iC s {a 1 ,x i ) 

-iC s (a, - C§ (a, yi)) 
ifj 



2 (d s C s (a,-)) (y) = 0, 



and similarly 



Pi[G 2 (x) - G 2 (y)} = ^(-C s (a,x 1 )+iC s (a,x_ i ) 



If £ = rj 3 R, we have that x — y + and 



+C S (a, ?/_,) - iC s (a, yi)) . 
2 



P/ [Gi (as) - Gi (»)] = -^(CjK^l-CiK^ 

+iC*«5 (a, y_i) - C 5 (a, yi)) 

3 

= --^(^C 5 (a,-))(.x)=0. 



and similarly 

3 

P e [G 2 (x)-G 2 (y)] = ^(-C^csiJ-tCaCca:-*) 

+C 5 (a, + iC,5 (a, yi)) 

= ^(a,C f (o,.))(y)=o. 

This concludes the proof of the lemma. □ 



We now turn to the singularity of fc s (Proposition 16 1: 

Proposition 33. Near the midpoint of a horizontal edge a £ Vc 6 , for x € {±1, H}, 
set a x := a + ^ € V c ju and oy e z := ( a 7 a x) G £c M - ^ erl ^ e function fc s (a, •) 
satisfies the relations 

p ^(e 1+i ) [/cj (a, a)] = P% 1+< ) [/c 5 (a, ai+i)] . 

P^( ei _i) [/Cj ( a , a )] = P^(ei_i) [/c 5 (a, , 

P|(e_ 1+1 ) [/c 6 (a, «) - 1] = p %_ 1+4 ) [/c 5 (a, a-i+i)] , 

P^(e_x_i) [/c 5 (a, a) - 1] = P^e-i-i) [/c 5 (a, a~i-i)] ■ 
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Proof. Set c := cos (|) and s := s in yy a nd rj := e" . Exact values of the coupling 
function Co that can be found in [KcnOO (see Figure 6 there) 



1 

4' 
i 

4' 



—i 


(i 






1 


i 




~ 4' 


i 


1 




7T 


i 


1 


4 " 


7T 



Co (0, 1) = -Co (0, -1) 
Co (CM) = -C o (0,-i) 
C (0, 2 + i) = C (0, -2 + i) 
Co (0, 1 + 2%) = C (1 - 2i) 
C (0, 2 - i) = Co (0, -2 - 2i) 

Co(0,-l-2i)=C (0,-l + 2i) 

Using these values and the definition of fc s , a straightforward computation gives 

'2 l + i N 



/c s ( a ; a l+i) 
/cj (o,ai_i) 

/Cj (OjO-l+i) 

/c 4 (a,a_i_i) 



77 c 



r? ^c 
jj f c 



7T 2 



2i 1 + i 
2i + 2 1 + i 



2 + 2i 1 + i 



V c 



7T 2 

2i 1 + i 
7 2~~ 



is 



2 

1+i 



/s 



2 1 + i 
n 2~~ 



If we compute the projections of these values on the lines associated with the medial 
edges e x , a straightforward computation gives 



P v 3 R [fc s (o,ai+i)] 
P»j3 R [f Cs (a, ai-i)] 
P »)R [fc s (a, a_ i +i )] 

Pr,R [fc s ( a ; a -l-i)] 



V — = P-n 
V2 " 1 



71 



2 + 72 
4 

2_+_v/2 
4 

2 + \/2 

~ 4 

2 + ^2 



which is the desired result. 



□ 



Let us now recall the result of Kenyon concerning the convergence of the function 



Co 



Theorem 34 (Theorem 11 in [KenOOj). As \z\ — > oo, we have 

f 3Re ( — ) +0 (tV) z = 2m+(2n + l)i:m,neZ, 
C o (0,z) = { VW / 

I iShn(i) +0 f j^WJ z = (2m + 1) + 2m : m,n G Z 
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From there, we can prove Theorem |25| 



Proof of Theorem \25\ By rescaling the lattice of the theorem above, one readily 
deduces that 



C I ^ I " 



id 



i 



0. 



2tt (z — a) 

uniformly on the sets {(a, z) : \a — z\ > e}. The proof of the theorem follows then 
from the definition of fc s ■ □ 
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